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Abstract 

Intracellular signaling often arises from ligand-induced oligomerization of cell 
surface receptors. This oligomerization or clustering process is fundamentally 
a cooperative behavior between near-neighbor receptor molecules; the proper- 
ties of this cooperative process clearly affects the signal transduction. Recent 
investigations have revealed the molecular basis of receptor-receptor interac- 
tions, but a simple theoretical framework for using this data to predict cluster 
formation has been lacking. Here, we propose a simple, coarse-grained, phe- 
nomenological model for ligand-modulated receptor interactions and discuss 
its equilibrium properties via mean-field theory. The existence of a first-order 
transition for this model has immediate implications regarding the robustness 
of the cellular signaling response. 
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I. INTRODUCTION 



Cell growth, differentiation, migration, and apoptosis are in part regulated by extracel- 
lular polypeptide growth factors or cytokines (Heldin, 1995; Stuart and Jones, 1995). As 
these molecules are unable to pass through the hydrophobic cell membrane, they have to 
bind to the extracellular domains of specific surface receptors in order to exert their effects. 
Much effort has gone into investigating the fundamental question of how the ligand-receptor 
interaction can trigger the proper intracellular signals. One popular hypothesis is that 
ligand-induced "clustering" of ligand-receptor complexes can be a key element in the proper 
activation of downstream signals. (Ashkenazi and Dixit, 1998; Bray and Levin, 1998; Heldin, 
1995; Germain, 1997; Lemmon and Schlessinger, 1994, 1998; Reich et a/., 1997; Sakihama 
et al, 1995). 

As an example of this line of reasoning, we consider the signaling cascade mediated by 
the binding of tumor necrosis factor (TNF) to the receptor TNF-Rl. Internally, the cyto- 
plasmic domain of TNF-Rl is "sensed" by a variety of adaptor proteins, namely TRADD, 
FADD, TRAF2, and RIP; this sensing leads eventually to NF-kB/JNK/SAPK activation 
and apoptosis. To accomplish the downstream signaling, an oligomerization of these adap- 
tor proteins is required (Ashkenazi and Dixit, 1998). One way to facilitate oligomerization 
is via construction of a molecular scaffolding via TNF-induced TNF-Rl clustering. It is 
known that TNF-Rl will not aggregate in the absence of TNF; this is due to the association 
of an inhibitor, "silencer of death domain" (SODD), which normally attaches to TNF-Rl 
cj^oplasmic domains and prevents receptor aggregation (Jiang et ai, 1999), or alternatively 
due to the receptor extracellular domains since spontaneous association of TNF-Rl has been 
observed in cells that express truncated receptors (Boldin et al, 1995; Vandevoorde et al., 
1997). TNF treatment, however, can bring two or more receptors into proximity via its 
multiple binding capacity (Jones et ai, 1990, 1992)]. This "proximity" might "squeeze" out 
SODD (Jiang et al., 1999), expose the cytoplasmic "death" domains to adaptor proteins, 
and thereby stabihze receptor clusters. Thus, a molecular scaffold/nuclei is generated to 
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initiate signaling. 

Over a longer time scale, the signaling messages can provide feedback to modify the 
capability of surface receptor clustering (Humphries, 1996; Wyszynski et al., 1997). This 
leads to a complex dynamical process involving both the intracellular signaling cascades 
as well as the surface receptor clustering. The self-organization made possible by these 
feedbacks has been intensively discussed for signaling cascades (see, e.g., Jafri and Keizer, 
1995; Barkai and Leibler, 1997). Much less is understood, however, regarding the role of 
receptor clustering. It is clear, though, that given the hypothesis that cellular signahng relies 
on the formation of receptor clusters, the temporal and spatial characteristics of clustering 
would certainly affect the process of signaling transduction. Thus, modeling the physical 
properties of receptor clustering is as important as modeling signaling cascades. 

Since clustering is due to an interaction between nearest-neighbor receptors, it is obvi- 
ously a cooperative process. Prom a physics perspective a system with this type of coop- 
erativity can exhibit a first order phase transition, corresponding to a jump in the surface 
density of ligand-receptor complexes. In the coexistence region of this transition, the surface 
will spontaneously segregate into two phases, dilute and dense. This first order phase tran- 
sition endows the signal transduction process with the ability to produce a digital signal in 
an analog world; this is independent of the details of intracellular cascades, instead, arising 
from the intrinsic cooperativity in ligand-receptor interaction. This has not been adequately 
addressed in the few models studied to date (Goldstein and Wiegel, 1983; Goldstein and 
Perelson, 1984; Riley et al, 1995; Coutsias et al, 1997; Shea et al, 1997). 

The purpose of this work is to introduce a phenomenological model for the TNF-TNFRl 
system to describe the onset of receptor clustering (phase separation). Specifically, we 
assume that clustering can be described by the statistical mechanics of a simple lattice 
Hamiltonian, incorporating the fundamental mechanism of a multimeric binding capacity 
for the hgand. We will calculate (via mean-field theory) a phase diagram and show that 
clustering will be thermodynamically favored for some range of ligand and receptor densities. 
Finally, we will do a simple Monte Carlo simulation of this system, showing that receptor 
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diffusion will lead rapidly to cluster formation in the relevant parameter range. We neglect 
the possibility that there exist long-time feedbacks to modify the clustering capacity, and 
we ignore some inessential details of the receptor-ligand interaction. More detailed models 
including these effects, as well as applications to other signaling systems, will be presented 
in the future. 



II. THE LATTICE HAMILTONIAN 

In our model, we treat the cell surface as a lattice with a spacing of the order of a 
few nm; this is the closest that neighboring receptors can get to each other. Each lattice 
site i has either one or zero receptor molecules, denoted as = 1 or 0. Our receptor has 
only two states: liganded or unliganded, and the interaction between receptor molecules 
is determined by their states. This "two-state" model is over-simphfied, yet we will see 
that it gives reasonable predictions for the phase diagram. A "state" label, = 1 or 2, to 
represent unliganded or liganded, then, can be assigned to each occupied receptor. We will 
further assume that the only ligands on the surface are those bound to receptors. If we let 
the chemical potential of the ligand be /^l and that of the receptor be hr, we then get a 
contribution to the effective Hamiltonian of the system 

H,{{n,t})^-Y^fi{ti)ni (1) 

i 

where /^.(l) = /ir and /x(2) = + + and —gi, is the binding energy between ligand 
and receptor. 

We should clarify the relationship between the parameters used here and those in real 
experiments. Using standard ideas (Changeux et ai, 1967), we notice that with only this 
term, the partition function can be factorized and reduced to a single site problem. 



fc=i 



(2) 



From this, we can immediately obtain the expectation values of the TNF-Rl concentra- 
tion in the liganded and unliganded states. These are assumed to correspond to the 



equilibrium condition of the following reaction (Corti et ai, 1994; Grell et al, 1998): 
TNF-R1'^™^+TNF ^ TNF • TNF-Rl'^™^ with a corresponding equilibrium dissociation con- 
stant: }™!;~il^|r'l:^i™^' = Kjr^ ^ 0.59 nM, where the notation "TNF-Rl^"^)" means TNF- 
Rl molecule distributed on the artificial membrane, and where the brackets [. . .]eq indicates 
the equilibrium concentration of the respective molecule. From this, we have e^^^^'^^'^^ = 
. To obtain the parameters individually, we might employ an "ideal gas law" for the 
hgand. This yields e^^^ = [TNF] ( .^^^^J^,^^ )'^', and = ksTln 
GOkBT, with h as Planck constant and mmf as the mass of TNF. 

We next add a receptor-receptor interaction term. This takes the general form 

H2{{n, t}) = Jija{ti, tj)ninj (3) 

Here, Jij = 1 only when i,j are nearest-neighbors and is otherwise (fig.|I]). The function 
a(ti,tj) indicates a "state" -dependent interaction energy between nearest-neighbor recep- 
tors, namely, a(l, 1) is the energy between two unliganded receptors, a(l,2) = a(2, 1) is 
the energy between one liganded and one unliganded receptors, and a(2, 2) is the energy 
between two liganded receptors. We note that in general, higher order terms might exist, 
especially considering the "trimeric" nature of the TNF ligand in our model problem. We 
have similarly neglected the details of the interactions of the ctyoplasmic domains, as per 
our earlier discussion. Our goal is to elucidate the basic idea regarding clustering in the 
simplest possible model, assured that adding more details will not change the basic notion 
that there exists a first-order transition due to the cooperativity. 

a(l, 2) and a(2, 1) are the interaction energies, for which we will use an effective binding 
strength qe of the order of Ql/^O, arising via one or two hydrogen bonds between receptors. 
It is important to realize that our simplified model does not treat explicitly the formation of 
multimers via the multimeric binding. Instead, it arbitrarily assigns the one ligand (binding 
two receptors into a dimer, e.g.) to one of the receptors and describes the dimeric binding 
as an attraction between a bound and an unbound receptor. Because of this, the model 
cannot distinguish between this relatively strong interaction and the subsequent much weaker 
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interaction between the dimers. In future work, we will show that this complication does 
not alter the basic picture presented here. 

As discussed above, in the TNF system there is probably a short-range and non-specific 
"excluding" interaction between two unliganded or two liganded (with different ligand 
molecules) receptors. For the sake of simplicity, we will assume that the repulsive energy is 
of the same order of the associative one, i.e., a(l, 1) ~ a(2,2) ^ —Qe- This assumption is 
not necessary, yet it greatly simplifies the mathematical task for analysis. 

The symmetry of a{ti,tj) allows us to introduce a simple matrix notation for the total 



Hamiltonian H1 + H2. If we use two-component vectors for the state-labeling: Tj 



for 



ti — 1, and Ti 



for ti = 2, then the Hamiltonian can be rewritten as 



H{{ni,Ti}) = /i(2)]Ti - Jy-nin^Tj^ 

i ij 



-1 1 
1 -1 



(4) 



Here [/^(l), A*(2)] is an 1 x 2 matrix. The simplicity of using this form of the matrix a{ti, tj) 



can immediately be seen if we make a transformation Ti — ^ 



1 - (Ti 



, with (Ti = ±1. Then 
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H{{ni, ai}) ^ -xY^rii- niGi + -qe Y1 JijUiUjaiaj 

i i ij 



(5) 



where x — -|- /x(2)]/2 is the "averaged" receptor chemical potential, and y — — 



/i(2)]/2 is directly related to the hgand concentration, e — \ r^n) ■ The partition function 

V d 

then reads 



exp[-/3i/({n,a})] 

{ni,cJi} 



(6) 



where Y,{rn;ai} means ensemble summation over the three different configurations {rij = 
0; rii — l,ai — ±1} on each lattice site, (3 — l/ksT with ks as Boltzmann factor and T as 
temperature. 
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If we define a new notation Ui — riiai, our model would be very similar to a spin-1 
antiferromagnetic (AFM) BEG model (Blume et al, 1971), 

Hi{ui}) ^ -xY.u'^i -yJ2ui + \gE JijUiUj (7) 

i i ij 

The origin of this AFM behavior is the "negative cooperation" between nearest-neighbor 
receptors, as we have imposed that a "proximity" of two unliganded or two liganded recep- 
tors will cost energy. Similar behavior might occur in the erythropoietin receptor (EPO-R) 
and the human Growth hormone receptor (hGH-R) systems (Heldin, 1995). This negative 
cooperation will give rise to an absence of clustering in extreme high/low ligand concentra- 
tion (i.e., y — > ±cxd), and thereby result in a "bell" -shape or window- like signaling response 
(Elliott et al, 1996). 

We should point out that this negative cooperation is not universal. In the case of EGF- 
R (epidermal growth factor receptor) system, a ferromagnetic (FM) behavior ("positive 
cooperation") is more likely, since there clustering requires two or more liganded receptors 
(Lemmon et al, 1997). Thus the higher the ligand concentration, the more the EGF-R 
cluster can be formed, and the EGF-EGFR signaling response behaves in a sigmoidal rather 
than a window-like pattern. It is clear that in both EGF-R and TNF-R (and hGH-R, 
EPO-R) systems, the ligand multiple binding capacity is the essential ingredient to induce 
clustering (of course one should consider the effect from receptor cytoplasmic domain as 
well). Which kind of cooperation (negative or positive) one should one consider depends 
on the details of the receptor-receptor interaction (also including the chemical modifications 
on receptor cytoplasmic domains), and needs to be be established experimentally. But, 
the essential feature of a first-order-transition-like behavior in receptor clustering is not 
dependent on the sign of this additional cooperativity. 

III. NUMERICAL SIMULATION 

To see if our model can generate clustering, we perform a Monte Carlo simulation on 
a square lattice with the standard Metropolis scheme. For simplicity, we fix the number 
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liganded and unliganded receptors and do not allowed these to fluctuate. Given the rather 
strong binding, this is not an important constraint. Furthermore, we allow motion only for 
individual recpetors and do not explicitly allow a clusters to move as a whole; this might 
not be the case in reality. The "jumping" probability for a receptor to move to another 
lattice site is determined by the Hamiltonian and obeys the detailed balance law. In detail, 
we pick a receptor at random and try to move it in a randomly chosen direction. The move 
is accepted if it lowers the energy and the move is accepted with probability e~^^^ if the 
energy increases. 

From fig.^, we immediately see that for a given receptor density, changing the ligand 
concentration moves the system from a non-clustering to a clustering phase. In this figure, 
the open, filled circles indicate liganded, unliganded receptor molecules, respectively. Note 
that the open and filled circles are arranged in an alternative way to form the cluster (i.e., 
inside a cluster, the nearest neighbors of the open circles must be filled circles, and vice 
versa). This implies that the equilibrium state (which must be translationally invariant), 
can be described via dividing the system into two interleaved sub-lattice systems: one sub- 
lattice is occupied by one species of receptor molecule (liganded or unliganded), and all its 
nearest neighbors belong to the alternative sub-lattice which is occupied by another species. 

To obtain more insight into the conditions where receptor clustering can take place, we 
next analyze the partition function via the mean-field approximation. 

IV. MEAN FIELD APPROXIMATION 

To proceed, we decouple the quadratic term in the Hamiltonian by introducing an aux- 
iliary Gaussian field and employing the standard Hubbard-Stratonovich/ Gaussian transfor- 
mation (see, e.g., Amit, 1993; Parisi, 1988), eqn.( [A3| ). The benefit of this transformation is 
to decouple the quadratic terms into linear terms such that we can sum over the ensemble 
configuration {{rii, cTj}) at each lattice site i independently. This yields (see Appendix for 
details) 
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Z = C V(f)e 



(8) 



with Si{{(j)}) = In l + 2z cosh{p[gE JiAi + v]) 



j 



where V(j) = Ylid(j)i, and C is a normahzation constant which does not affect the thermo- 
dynamic properties of the partition function. The new field ranges from — oo to +00, 
y = — (/iL + Ql)/'^-, and z = e^^ = e^^^^^^^^'^^^'^^'^'^ is related to the receptor chemical poten- 
tial which remains to be determined (in terms of the receptor density). The first term in 
eqn.(D is related to the interaction energy between nearest-neighbor lattice points, whereas 
the second term is related to the entropy arising due to the available configurations on an 
individual lattice site. 

In mean-field theory, we try to determine a "homogeneous" saddle point approximation 
for the partition function. For our system, the negative cooperation (i.e. the AFM nature) 
suggests that the system might prefer having neighboring sites in oppositely liganded states. 
Thus, we separate the lattice into two interleaved sub-lattice systems: all nearest-neighbors 
of a lattice site belong to the alternate sub-lattice (fig.|l]). We then assign two "uniform" order 
parameters, 0± to each sublattice. After this assumption, the exponent of the Boltzmann 
factor in the partition function [eqn.(|)] now becomes — > ^ [l3gED(j)^(j)^ + S'(0+, 0-)], where 
N is the number of total lattice sites, S'(0+, 0„) = Y.k=± ln[l + 22; cosh(/3[(7£;D0fc + ?/])], and D 
is the number of nearest neighbors, which depends on the structure of lattice. For instance, 
a square lattice yields D = 4, whereas a honeycomb lattice yields D = ?>. 

Next, we minimize the free energy by varying (f)±. The variation yields the "saddle point" 
equation 








. Working this out explicitly, we find a self-consistent equation for (j)± 



(9) 



l + 2zcosh.[(3{gED(t)^ + y)] 



with the free energy density 
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z) = -^0+0- - ^ E + 2zcosh{(3[gED^k + 2/])] (10) 

k=± 

Finally, the mean-field receptor density is given by (n) = —df{(j)+,(j)-,z)/dfiR. Explicitly, 
we have 

^ j2 zcosh{(3[gED(i)k + y]) ^^^^ 
1 + 2zcosh{p[gED(j)k + y]) 
We can therefore determine the receptor chemical potential, x (or equivalently z) in terms 

of (n). Thereafter, we can rewrite the free energy density in terms of (n), 0+, and 0_. 



V. THE ONSET OF CLUSTERING 

There is no closed form solution for eqn.(^. To get some analytical information, we 
define 0± = m ± e and, with U{w) = ifarctSSS^V '^^^^ 

oo 2k 

^=-j:^,u^''\^) (12) 

-EpI^T^f'^'-'M (13) 

where U^^\w) = d^U{w)/dw^. The basic idea of separating out the e dependence is that 
solutions with non-zero values of e represent phases in which the proximity of neighboring 
receptors gives rise to alternating ligand binding. For very small receptor densities, there are 
few neighboring receptors and hence we expect to find a unique solution of the mean-field 
equations with e = 0. In fact, it is clear from eqn. ([T3|) that there is a solution with e = for 
all values of the parameters, but at larger densities, there may be other more stable phases. 
The goal of our analysis will be to understand the general structure of the phase diagram 
and then to obtain more quantitative detail by numerical means. 

To proceed, let us assume that e is small and solve eqns.(|l^), (|13]) to order e^; 



m = mo — rriie (14) 



with mo = -f/(°H 



mo J 



= 2![l + f/a)(mo)] ^''^ 
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mif/(2)(mo) - |f/(3)(mo) ^ ^ 

Using the relationship given above for (n), It is easy to verify that U'^^\mo) = l3gED\{n) — 
mg], f/(2)(mo) = -mo{(3gEDf[l - 3(n) + 2m2], and 

f/(3)(mo) = ([l - ^n) + 6m2] [(n) - m^] - Sm^ [1 - (n)]) 



We must consider separately the cases where the denominator of eqn. (|Tq) is positive 
or negative. Let us first imagine it is positive, Then, the existence of a non-trivial solution 
of eqn. (|l6|) requires that {(3gED[{n) — wiq] — 1} > 0. At small {n) this condition will 
clearly fail and we will have only the trivial solution. Also, this condition will fail at (n) 
close to 1 for large enough \y\. We can see this by comparing the equation for itlq with the 
expression for (n). Note that if y is large enough such that the hyperbolic functions can 
be replaced by exponentials, we have |mo| = (n), and the above expression can be replaced 
by {(3gED[{n) — {n)"^] — 1}; this is negative for the stated condition. As we cross a line in 
parameter space such that this factor changes sign to positive, there will be new solutions 
at non-zero and the one at e = becomes a local maximum of the free energy. This 
emergence of a double-well structure, with a continuous growth the non-zero solution 
indicates that the system exhibits a second-order phase transition. 

We must next take into account the possibility that |mif/(^^(mo) — |?7*^^-*(mo)| < 0. 
Having the denominator cross zero gives rise in our current approximation to a large value 
of e which thus invalidates the neglect of higher-order terms. Typically, the higher-order 
terms will stabilize the system at some finite value of e, which thus appears "spontaneously" 
as some parameter threshold is crossed. This is a first-order phase transition, or equivalently 
a triple- well structure for the free energy. If the local minima (for zero and non-zero e^) have 
equally low free energy density, the system can exist in a mixture of the two phases. As 
we will see, the two coexisting phases differ in their receptor density. Finally, the points 
where both {l3gED[(n) — mg — 1] = 0} and |mi[/'^^)(mo) — |f/'^^^(mo)| = 0, are "critical 
end-points" points, since they correspond to places where a second-order transition line 
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ends at a first-order line. A diagram of this beliavior, generated by tlie numerical solution 
of the mean-field equations, is given in fig. ^. 

For a given ligand concentration, we can find the phase coexistence lines arising due to 
the first-order phase transition. This is done by finding two solutions (solved with differing 
values of e^) of the mean field equations and then fixing z (as a function of y) by requiring 
that they have equal free-energy 

/(0W,0('^),^) = /(0L^0j\z) (17) 

where (\)~^ ^ are the order parameters for the dense condensed phase and (\)^'^'^ the (equal) 
ones for the dilute phase. For the condensed phase, the receptor density is close to unity for 
reasonable values of the cooperativity parameter (SqeD. The workings of this system as far 
as signaling is concerned is shown in fig. ^ Assume there is some fixed value of the receptor 
density. As the ligand concentration is increased, we will cross the phase transition boundary 
and the receptors will segregate into a condensed phase and a dilute one, corresponding to 
the two co-existing mean-field solutions. Under our basic hypothesis that signaling is effected 
by having dense clusters, the response will exhibit a sharp jump at a specific threshold ligand 
concentration. Similarly, as the ligand concentration becomes too high we cross back to the 
uniform receptor density state and signaling ceases. That is, we have a ligand concentration 
"window" for receptor clustering. 

As can seen from the figure, the "clustering" window will cease to exist below some 
minimal receptor density, as we never enter the phase coexistence region. By symmetry, this 
minimal density can be found by solving the mean field equations for y = where m = 0. 
This leads after some algebra to the self- consistent equations 

, „ (.PqeD^V^ _ 1 

(^)Lt = ni , (18) 

cosn\i3gEDe\ — 1 

with e = (n)^j-„ tanh[/5(y'£;-De] 

The numerical solution of these equations is presented in fig.^ As the cooperativity pa- 
rameter is increased, the minimum density which will support a clustering window goes 
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rapidly to zero. For the TNF-TNFRl cluster, it has been speculated that the structure 
of cluster is a honeycomb-like lattice (Bazzoni and Beutler, 1995; Naismith et al, 1995, 
1996), which implies the number of nearest-neighbor D = 3. If we use our rough estimate 
Qe — 5'l/IO ~ 6kT, we find that ^ lO'^/a^. Here qq is length scale of the lattice 

spacing. If we take Oq ~ 1 nm, on a cell with surface area 100/xm^, this estimate yields 
a requirement for less than 10^ TNF-Rl molecule distributed on the cell surface. Given 
that an average number of expressed TNF-Rl on cell surface is ~ 2000, we find that the 
cell operates within the desired part of the phase diagram and hence should exhibit strong 
sensitivity to the application of TNF. However, we should point out that this estimate is 
very rough, as we have made a number of simplifying assumptions and this issue needs to 
be re- visited with a more precise model of the receptor interactions. 



VI. DISCUSSION 

We have presented a simple model for signal transduction via receptor clustering, based 
loosely on the TNF-TNFRl system. Our basic idea is simple. The interaction between 
receptors can lead to a first-order phase transition with a discontinuous jump in the receptor 
density as a function of the receptor chemical potential and/or the ligand concentration. 
Turning this around, this implies that the receptor system will spontaneously phase separate 
for a range of ligand concentrations. This fact about the thermodynamic equilibrium state 
will lead under reasonable kinetic assumptions to the rapid formation of receptor clusters. 
Assuming that these clusters are necessary for the signal to proceed downstream has the 
immediate consequence that the system exhibits a strong robust response independent of 
any details of the intracellular signaling cascade. This might provide a simple solution to 
the problem faced by biological evolution of how to get digital response in an analog world. 

From a physics perspective, there is nothing very surprising about our phase diagram 
findings. The idea of a "lattice" Hamiltonian with intrinsic "cooperativity" has been pro- 
posed before (Changeux et al, 1967), and on general grounds models of this sort can be 
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expected to have first order phase transitions. What is new here is the connection of the tran- 
sition to signahng via the idea of receptor clustering. This connects nicely with increasing 
evidence that clustering is "universal" among many types of receptor classes. 

In our model, we have ignored more-than-two receptor interaction, and relevant internal 
chemical degrees of freedom (such as the dissociation of SODD in the TNF-Rl system). 
We do not expect these detailed considerations to change the overall picture, but a more 
sophisticated model will be needed to make more quantitative estimates of ligand thresholds, 
cluster structures and most interestingly, clustering dynamics. We hope to report on these 
issues in the future, as well as on the extension of our models to other ligand-receptor 
systems. 

Finally, it would be important to extend our work to later-stage dynamics, as that would 
allow the consideration of processes such as adaptor protein-mediated receptor internaliza- 
tion, cytoskeleton-assisted cluster stabilization, receptor affinity regulation, receptor cross 
talk, and adaptation (Barkai and Leibler, 1997; Hahn et al, 1993; Humphries et al, 1996; 
Holsinger et al, 1998; Luo and Lodish, 1997; Stewart et ai, 1998; Sundberg and Rubin, 
1996; Valitutti et al, 1995; Wyszynski et al, 1997). Other possible extensions might in- 
volve the inclusion of spatial fluctuations, the explicit treatment of external perturbations 
(Shoyab and Todaro, 1981), the local heterogeneity of the micro-environment (Bean et al, 
1988; Ward and Hammer, 1992), or fluctuations of hgand concentration; all of these issues 
have been neglected here. 
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APPENDIX A: THE GAUSSIAN TRANSFORMATION 



The identity 



dxexp 



1 2 • 

X + ISX 

^9 



C X exp 



2 



(Al) 



with g > 0, i = and C as some constant, can be generahzed to 



Y[ dxi exp 



jk j k 



= C X exp 









(A2) 



as long as J is a symmetric positive definite matrix. Thus, eqn.(H) can be obtained by 
Z = expl(3[xY^ni + yJ2niai~^Y^Jijniainjaj 

{ni,(Ti} i i ij 



(A3) 



^0 /n<^ie 



_{nj,(7j} 



= / n '^j-e^ X n 1 1 + 2^ cosh jjk<i>k + 

•' 3 i I Ik 

where i = 0^ = z6'fc, = In [1 + 2^; cosh(/5[5f£; Efc -'jfc^fc - 

integral constants. 



and Cg, are 
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FIG. 2. Monte Carlo simulation with Metropolis scheme. Here we test the model under a fixed 
receptor density but different ligand concentrations. In both upper and lower panels, the left figure 
represent initial conditions and the right are results after 10^ Monte Carlo steps. The open, filled 
circles indicate liganded, unliganded receptor molecules, respectively. There is no stable cluster 
formation in the upper snapshot, whereas the clustering is stable in the lower one. Here we use 
gE = GksT, density of liganded receptor: upper plane, 0.001, lower plane, 0.03, and density of 
unliganded receptor: upper plane, 0.059, lower plane, 0.03. 
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FIG. 3. Numerically computed phase diagram, showing that there is a pair of second-order 
lines each of which ends on the first-order transition curve. (C=critical point; CE=critical end- 
point). The phase to the lower-right has e 7^ 0. Here we used qe = QksT and D = 3. To show the 
symmetry, we plot the ligand concentration in logarithm unit, normalized with respect to the dis- 
sociation constant K^J^^ = 0.59 nM. Here the chemical potential is related to receptor density. 
In fig. ^, we convert the receptor chemical potential into the molecular density. 
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FIG. 4. The phase diagram now shown as a function of the receptor density (n) related to 
the hgand concentration. To show the symmetry, we plot the ligand concentration in logarithm 
unit, normalized with respect to the dissociation constant K^J^^ = 0.59 nM. The region inside the 
solid lines is the co-existence region where states of high and low density co-exist. As the ligand 
concentration is altered so as to cross one of these lines, the receptors will spontaneously cluster 
and thereby allow signaling to occur, (n)^]^ is the minimal receptor density for clustering. For 
this set of parameters, clustering will occur even for very small overall receptor density. Here oq is 
the length scale for the lattice spacing. For surface receptor molecule such as TNF-Rl, we might 
take ao ~ 1 rim. 
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FIG. 5. The variation of as a function of the association energy qe- We find that 

('^)mL rapidly approaches zero once gsD > ISksT. If we assign £) = 3 ~ 4, this energy scale 
corresponds to a single hydrogen bond. Here oq is the length scale for the lattice spacing. For a 
surface receptor molecule such as TNF-Rl, we might take oq ~ 1 nm. 
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